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Abstract 

The gluon dipole operator, On, has received much attention recently be- 
cause it can have a large coefficient in some SUSY extensions of the standard 
model. We find that the commonly used matrix element of On of Bertolini, 
Eeg and Fabbrichesi is in rough (but accidental) numerical agreement with an 
estimate based on dimensional analysis. 



One of the operators in the low energy effective Hamiltonian responsible for \AS\ = 
1 weak decays is the gluon dipole operator, which can be written as (-f^e// = CiiOu), 

Heff = Cu^s[m,R + m,L\T^G^^a^, + h.c, (1) 

where C^" is the gluon field strength tensor, and L, i? = (1 =1= 75)/2. In the standard 
model, the coefficient Cu is sufficiently small to make the effects of this operator 
negligible. Beyond the standard model, however, this operator can have a large 
coefficient and it becomes important to estimate its matrix element. The operator 
has received much attention recently, in connection with e'/e 0. 

The importance of this operator for kaon decays, and in particular for the analysis 
of e'/e dates back to the Weinberg model of CP violation in the early 80's. It was 
shown back then, that one had to pay particular attention to the chiral properties of 
the operator [Q in order to obtain physical amplitudes that obeyed the FKW theorem 
i- 

The gluon-dipole operator transforms as (3, 3) or (3, 3) under chiral rotations and, 
in the standard model, it is proportional to the light quark masses. It is well known 
that there are no operators with these transformation properties in the lowest order, 
0{p^), weak chiral Lagrangian |^, §[. At next to leading order, 0{p'^), there are 
several operators with the desired properties. 

The strong interactions of pions and kaons are described to order 0{p'^) in chiral 
perturbation theory by the Lagrangian of Gasser and Leutwyler The ingredients to 
construct this Lagrangian are the non-linear function S = exp(2i0//) which contains 
the octet of pseudo-Goldstone bosons and that transforms as S — > KEL'^ under 
SUl{3) X SUfi{3). Interactions that respect chiral symmetry are constructed in terms 
of derivatives of S. Explicit chiral symmetry breaking due to the non-zero light quark 
mass matrix M = diag{mu, rud, m^) is introduced through the factors S and P, 

s = x^s + stx 

P = z(x1"S-Stx) (2) 

For our present purpose it suffices to take x = 2BqM, where the parameter Bq is 
proportional to the quark condensate < qq >, and relates the current quark masses 
to the meson masses: 

< >= -fn^o, m\ = Boinis + rrid), ml = Bo{mu + rrid) (3) 

The leading order weak chiral Lagrangian transforming as {8l,^r) under chiral 
symmetry was first written down by Cronin 0] and has only one term (We use the 
notation = iE'^'D^E), 

£SJ) = C2Tr(A6L2) (4) 

To introduce explicit chiral symmetry breaking due to the non-zero light quark mass 
matrix M into the chiral Lagrangian one pretends that the mass matrix transforms 
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as M — s> RML'^ and constructs operators with the desired transformation properties. 
For the dominant (8^, 1r) weak operators it is well known that there is no mass term 
at order p^, in accordance with the FKW theorem. 

The operators that occur at next to leading order in the weak chiral Lagrangian 
have been written down in Ref. |0. From their general list, those that could corre- 
spond to chiral realizations of the operator On are the ones that contain one factor 
of the quark masses, they are (only five are independent) 

C^^ = ^ioTr(^A6{5,L2})+^nTr(^A6L^SL^)+Ei2Tr(^A6L^)Tr(^{L^5}^ 

+ Ei3Tr(^A6s)Tr(^L2) +Ei4Tr(^A6L2)Tr(^s) +Ei5Tr(^A6«[P,i:2]) (5) 

From the point of view of the transformation properties of On under chiral symmetry 



all these terms are equally valid. It has been argued in Ref. |TT| that the precise form 
of the short distance operator in Eq. |I| requires a chiral representation in which Ag 
and the quark mass matrix appear next to each other. This requirement reduces the 
possible operators in Eq. ||to the ones multiplying _Ei3 and the combination Eiq — Ei^. 
Using Eq. ^we find for the matrix element of On, 

2 /2 / 

On < 7r+7r-|On|A-° > = j^{m\ - ml)Uml + 2ml)Eio + mlEii 

J IT ^ 

,2 r,^2 \ 7-1 I I o»^2 N rp , ^2 



+ {Ami - 2m^)Ei3 + {mi + 2mj,)Eu + mi^E.^j (6) 

From this expression it is clear that the matrix element, in general, does not vanish 
in the limit 0. 

Although the framework of Eq. |^ is completely general, it does not tell us the 
size of the coefficients Ei. We can estimate the size of the Eiq, for example, that 
is needed to match the operator On by using naive dimensional analysis For 
this purpose we write down a Lagrangian with the minimal number of fields that is 
contained in the term proportional to Eiq (other terms are related to this one by soft 
pion theorems), C = gM^'xdfj^Kd^n and match to On following Weinberg 0, 
with the result 

I I ~ - 

V / NDA ' 

From this we write 



On < vr+vr |On|A'° > ~ (mi - m: 



2.m^ Oi 



NDA ' ^' U 167r2 47r 



To evaluate this expression one would use a value of Qs ~ v47r. Note that an equiva- 
lent expression is obtained if one uses the term proportional to i?i3 in Eq. instead. 
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Bag Model Estimate 



An explicit calculation of < 7r|0ii|i^ > within the MIT bag model [TO], sup 



plemented with a soft pion theorem, led to the estimate of Donoghue and Holstein 

Cii < TT TT |C»ii|K > = -Cn— — — — — (9) 

J MIT oZtt-^ 2/^ 

Eq. ^, is a trivial rescaling of the actual calculation in Ref. @. The factor Aj^n is 
obtained numerically from the Bag model, A^-k = 0.4 GeV'^, and the last factor in 
Eq. ^ introduces the suppression required by the FKW theorem. In accordance with 
power counting they choose A ~ 1 GeV, corresponding to a matching of On into an 
chiral Lagrangian such as the term that multiplies Eiq in Eq. ^. Numerically, 
we find that 

Cn < 7r+7r-|Oii|is:° ^ 0.8 (Cn < 7r+7r-|Oii|ir° (10) 

/ MIT V / NDA 



Chiral Quark Model Estimate 



More recently, Bertolini, Eeg and Fabbrichesi have used a chiral quark model 
supplemented with some matching conditions to estimate that ITTl , 



Cn < 7r+7r-|Oii|fs:° > 



CQM 



^2/ X 2 ^11 

— {ms-md)m^-— 



11 

< qq >G 



There are two points that we want to stress about this expression. First, the last 
factor in Eq. |lT| is the model dependent quantity that arises from their chiral quark 
model, and that we take at face value. Second, the factor of arises from the 
requirement that the short distance operator On match into an order p'^ weak chiral 



Lagrangian of the form [p 



£ ~ Tr 



(12) 



Comparing this with the general form, Eq. ^, we see that the only term that is retained 
is that proportional to the operator whose coefficient is Eiq — Ei^. As mentioned 
above, the requirement that the quark mass matrix M and Xq appear next to each 
other selects this term plus the one whose coefficient is £^13. This latter one is dropped 
in Ref. |Tl| because all products of two traces are suppressed in their model. 

however, is very similar to our dimen- 



The specific numerical result of Ref. ||TT 



sional analysis estimate. This happens because the m^/m^- suppression that occurs 
in the weak chiral Lagrangian operators that occur in their matching is compensated 
by the large numerical coefficient, 11, in Eq. [Tl|. Numerically, 



Cn < 7r+7r-|Oii|ir° > 



CQM 



1.4 [Cn < 7r+7r-|Oii|fs:° > 



NDA 



(13) 
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We see that all three estimates are numerically very similar and in agreement with 
each other within the uncertainty of each approach. However, it is clear that the 
numerical agreement with the chiral quark model result of Ref. [O] is accidental. 



Beyond the Standard Model 

In all the cases that we have discussed, we have considered the matrix element of 
the operator written as in Eq. |^. In this form it appears that the operator vanishes 
in the chiral limit being proportional to the light-quark masses. This was, in fact, an 
important ingredient in the matching to a corresponding chiral Lagrangian. However, 
in some models of interest, this is just an artifact of the normalization; the coefficient 
Cii goes as mj^. In this case we want to construct a low energy meson Lagrangian 
that transforms as (3, 3) or (3, 3) but that is not proportional to the light quark 
masses. 

It is possible to write a term without derivatives, 

C = Tr(j:^h + h^i:^ (14) 

that does not contribute to K ^ mi amplitudes once tadpoles are properly sub- 
tracted. The leading order amplitude arises from a Lagrangian with two derivatives. 
In terms of the matrix h whose only non-zero entry is /123 = 1 we can write, for 
example. 



C = ^^iTr(^/iD^S^D^SS^j + h.c. 

+ (7^2Trfe/i + /its)Trfz}^StD^s) (15) 



The first term in Eq. yields a matrix element proportional to m^, but the second 
term does not. We find, 



Cn < vr+7r-|0n|i^° >= 9N2^ml + 0{ml/ml) (16) 



Although we cannot compute 9^2-, we can again estimate it with naive dimensional 
analysis following Weinberg. Noting that the operator does not contain terms with 
only two fields, we find 

gN2-^m , (17) 



resulting in 



167r2 



< TT+TT-IOnli^^ > ) ~ m,m^^ (1^ 



This result is equivalent to Eq. ^, differing only by factors of order one that cannot 
be accounted for with naive dimensional analysis. 
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